SMOLUCHOWSKI'S COAGULATION EQUATION: UNIQUENESS, 
NON-UNIQUENESS AND A HYDRODYNAMIC LIMIT FOR THE 

STOCHASTIC COALESCENT 



JAMES R. NORRIS 

Abstract. Sufficient conditions are given for existence and uniqueness in Smolu- 
chowski's coagulation equation, for a wide class of coagulation kernels and initial 
mass distributions. An example of non-uniqueness is constructed. The stochastic 
coalescent is shown to converge weakly to the solution of Smoluchowski's equation. 



1. Introduction 

Coagulation of particles, in pairs and over time, within a large system of particles, is 
a phenomenon widely observed and widely postulated in scientific models. Examples 
arise in the study of aerosols, of phase separation in liquid mixtures, in polymerization 
and astronomy. Typically, it is argued that the rate at which pairs of particles 
coagulate depends, for physical reasons, in a given way, on some positive parameter 
associated to each particle, such as size or mass. In the models we shall consider, it 
is further argued that the effects of spatial fluctuations in the size or mass density 
are negligible — for example, by supposing that the particles perform independent 
random motions on a time scale faster than the process of coagulation. 

The first mathematical model of this sort of process was proposed by Smoluchowski 
S16|] in 1916, see also Chandrasekhar ||Cha43|| . Smoluchowski argued that particles 



V 



of radius r would perform independent Brownian motions of variance proportional 
to 1/r, so pairs of particles of radii ri and r2 would meet at a rate proportional to 

('^i + r2){l/ri + l/r2). 

Expressed in terms of masses, this leads to the coagulation kernel 

i^(x,y) = (xl/3 + yV3)(^-l/3^^-l/3) 

for particles of masses x and y. Then, making some implicit assumptions about er- 
godic averages, Smoluchowski wrote down the following infinite system of differential 
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equations for the evolution of densities fi{x) of particles of mass x = 1,2, 3, 



d 



^ x—l oo 
y=l y=l 



Here, the first sum on the right corresponds to coagulation of smaller particles to 
produce one of mass x, whereas the second sum corresponds to removal of particles 
of mass X as they in turn coagulate to produce larger particles. 

More generally, in other models, such systems of equations are considered for many 



different coagulation kernels K, see for example ||Ald|| . Also, analogous integro- 
differential equations are considered which allow for a continuum of masses x. It 
is known by now that, for a suitable initial mass distribution ^q, Smoluchowski's 
original equations have a unique solution. Much progress has been made in deter- 
mining when existence and uniqueness holds for more general coagulation kernels, see 



McL62 


1, 1 


[VlcL64 


1, 


Whi80|, 




Hei92 



Hei92|| for discrete mass distributions. Never- 



theless many fundamental questions remain open, even for certain coagulation kernels 
studied extensively in applied sciences. 

In this paper we give some new positive results on the existence and uniqueness 
problem for Smoluchowski-type equations. In particular: 

• we prove existence of solutions for continuous coagulation kernels K such that 

K{x,y)/xy^O as {x,y) ^ oo 

extending a result of Jeon ||Jeo|| for the discrete case; 

• we prove local existence and uniqueness of solutions when K{x,y) < ip{x){p{y) 
for some continuous sublinear function (p : E ^ [0, oo), provided that the initial 
mass distribution /xq satisfies 

(f{xYfj.o (dx) < oo; 

(0,oo) 

• this allows us to treat the case where K{x, y) blows up as x or y 0, also 
to prove uniqueness in some cases when the mass distribution has no second, or 
even first, moment; 

• we can do without any local regularity conditions on K; 

• we do not have to assume that the initial mass distribution is discrete, nor that 
it has a density with respect to Lebesgue measure. 

We also construct in §3 an example of a coagulation kernel K and an initial mass 
distribution hq, such that Smoluchowski's equation has at least two distinct solutions, 
both of which are conservative, in the sense that 



for all t. 



xjjt (dx) = / x/io (dx) < oo 

(0,oo) J(0,oo) 
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Then in §4 we consider a stochastic system of coagulating particles, where particles 
of masses x and y coagulate at a rate proportional to K{x, y). We show that, when we 
can estabhsh uniqueness in Smoluchowski's equation, the particle system, suitably 
normalized, converges weakly to the solution of the deterministic equation. Thus 
we obtain a statistical derivation of Smoluchowski's equation. This goes some way 
towards resolving Problem 10 in Aldous' survey article | |Ald| . 



2. Existence and uniqueness in Smoluchowski's coagulation equation 

Let E = (0, oo) and let K : E x E ^ [0, oo) be a symmetric measurable function, 
the coagulation kernel. Denote by M = M^; the space of signed Radon measures on 
E, that is to say, those signed measures having finite total variation on each compact 
subset of E. Denote by M"*" the set of (non-negative) measures in M. If /x G M"*' 
satisfies, for all compact sets B (1 E 

K{x,y)fi {dx)fi{dy) < oo, 

BxE 

then we define L{fi) G M by 

(/, L{f^)) = \ i {/(a; + y)- f{x) - f{y)}K{x, y)fi{dx)fi{dy) 
^ Jexe 

for all bounded measurable functions / of compact support. 

We consider the following weak form of Smoluchowski's coagulation equation 

(2.1) /^t = /io + / L{f^s) ds. 

Jo 

We admit as a local solution any map 

fit: [0, T) ^ M+ 

where T G (0, oo], such that: 

(i) we have 

/ xlx<i iJ,o{dx) < oo; 
Je 

(ii) for all compact sets B C E, the following map is measurable 

t^/it(5):[0,T)^[0,oo); 

(iii) we have, for alH < T and all compact sets B <Z E 

K{x,y)fis{dx)fis{dy)ds < oo; 



JBxE 
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(iv) for all bounded measurable functions / of compact support and also for f{x) = 
xlx<i, for alH < T 

(2.2) (/,/i*) = (/,/io)+ [\f,L{^.))ds. 

Jo 

In the case T = oo, we have a solution. 

The condition that, for f{x) = a;la;<i, we have (/, /Uq) < oc and that ( ^.21 ) holds is 
a boundary condition, expressing that no mass enters at 0. We obtain an equivalent 
condition on replacing / by any non-vanishing sublinear function E [0, oo) of 
bounded support, which is linear near 0. A function f : E ^ [0, oo) is sublinear if 

/(Ax) < A/(x) for all xeE,X> 1. 

Note that such a function / is always subadditive: 

fix + y) < f{x) + f{y) for all x,y e E. 

Hence (/, L{fi)) < for all G M+. Note also that, if ip : E ^ [0, oo) is any sublinear 
function and if 

{nx(f{n~^), < X < n^^, 
(p{x), < X < n, 

0, X > n, 

then (pn{x) I ^p{x) for all x, and ipn is sublinear of bounded support, linear near 0. 
So, for t < T 



Jo 

Hence, using monotone convergence on the left and Fatou's lemma on the right 

(2.3) (y?, /io) > (y?, /ii) - / {(fi, L{fXs)) ds. 

Jo 

In particular, {(f, fit) is non- increasing in t. In particular, the total mass density 

xfit (dx) 



is non-increasing in t; if it is finite and constant, we say that {nt)t<T is conservative. 
Throughout this section we make the basic assumption that 

(2.4) K{x,y) < !f{x)(f{y) for aA\ x,y e E 

where (p : E ^ {0, oo) is a continuous sublinear function. We also assume that the 
initial measure /iq satisfies 

(2.5) (V9,yUo)<00. 

We call any local solution {fit)t<T such that 



/ /is) ds < oo for alH < T 
Jo 
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a strong local solution. 

Here is a summary of what is known so far about existence, uniqueness and con- 
servation of mass in Smoluchowski's equation. The picture is more complete for 
discrete mass distributions — that is when /io is supported on N. Then, provided /iq 
has a finite second moment. Ball and Carr |BC90| proved existence and mass con- 



servation when K{x, y) < x + y, Heilman | |Hei92| | added uniqueness under the same 



hypotheses. Jeon [|Jeo|| has recently proved global existence when K{x,y)/xy — as 
{x,y) — >■ oo. McLeod |[McL62|| long ago proved local existence when K{x,y) < xy. 



For general mass distributions /iq, less is known. Dubovskii and Stewart ||DS96|| have 
shown existence and uniqueness provided /iq has an exponential moment and a con- 
tinuous density with respect to Lebesgue measure, and provided K is continuous 
with K{x, y) < I + X + y. Recently, Clark and Katsouros ||CK|| proved existence and 



uniqueness for a particular choice of kernel which blows up when x or y is small. 



Theorem 2.1. Assume conditions ( |2.4| ) and ( |2.5| ). If {fit)t<T CL^'d {i^t)t<T o,^^ local 



solutions, starting from /iq, and if {fJ,t)t<T is strong, then fit = Ut for all t < T. 
If (p{x) > ex for all x, for some e > 0, then any strong solution is conservative. 
Moreover, if {ip^,fio) < oo, then 

(i) there exists a unique maximal strong solution {fit)t<c(fio)^ with({fio) > (y9^,/io)~^, 

(ii) if ip'^ is suhlinear or if K{x,y) < (p{x) + ip{y) for all x,y E E, then C(/^o) = oo. 

The proof will occupy the remainder of this section. The method is to find an 
approximation to Smoluchowski's equation by a system depending on K and ip only 
through their values on a given compact set. The idea of the approximation may be 
readily understood in terms of the stochastic coalescent, for which a directly anal- 
ogous approximation is discussed in §4. Moreover, the finite particle interpretation 
explains certain crucial non-negativity statements, which are given less transparent 
analytical proofs below. 

We remark that this theorem provides examples where uniqueness holds, even when 
the solution fails to be conservative, in the trivial sense that the total mass density 
is infinite. We have not yet found an example of a strong solution which has finite 
initial mass density and which fails to conserve mass. The example of §3 shows, on 
the other hand, that uniqueness can fail while solutions remain conservative. 

Let B (1 E he compact. We will eventually pass to the limit B ^ E. Denote by 
Ms the space of finite signed measures supported on B. Note that ip is bounded on 
B. We define L^:MbxM— >M_BxMby the requirement 

((/,a),L^(/i,A)) = 

7T / {fix + y) ^x+y£B + dfix + y)lx+y^B ^ fix) - f{y)}K{x,y)n{dx)fi{dy) 
^ J ExE 

+ A / {aif{x) — f{x)}ip{x)fi{dx) 
Je 
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for all bounded measurable functions f on E and all a G M. Here we used the notation 
{{f,a),{fi,X)) = {f,fi) + aX. 
Consider the equation 

(2.6) (/it, At) = (/io, Ao) + / L^{^sAs)ds. 

Jo 

We admit local solution any continuous map 

t ^ {fit, Xt) : [0, T] ^ Mb X R 

where T G (0, oo), which satisfies the equation for all t G [0,T]. When [0,T] is 
replaced by [0, oo) we get the notion of solution. 

Proposition 2.2. Suppose fiQ G with /iq > and that Aq G [0, oo). The equation 
( |2.6| ) has a unique solution {fit,Xt)t>o starting from (/io,Ao). Moreover fit > and 
Xt>0 for all t. 

Proof. Our basic assumption ( p.4|) remains valid when ip is replaced by + 1, so we 
may assume without loss that f > 1. By a scaling argument we may assume, also 
without loss, that 

(V9,/Xo) + Ao < 1 

which implies that 

llyUoll + |Ao| < 1. 

We shall show, by a standard iterative scheme, that there is a constant T > 0, 
depending only on ip and B, and a unique local solution {fit, Xt)t<T starting from 
{fiQ, Ao). Then we shall show, moreover, that fit^^ for all t G [0,T]. 

First of all, let us see that this is enough to prove the proposition. If we put / = 
and a = 1 in (^.61), we obtain 



4:Xt = T: I ip{x + y)U+y^BK{x,y)fit{dx)fit{dy) + Xt / ip{xY fit{dx) . 
at Z .IexE 



So, since /xt > 0, we deduce Xt > for all t. Next, we put f = (fi and a = 1 to see 
that 

— ((y?, fit) + Xt) = - / {ip{x + y)- (p{x) - (p{y)}K{x, y)fit{dx)fit{dy) < 0. 
Hence 

WfirW + \Xt\ < {v,IJ't) + Xt < fi^o) + Aq < 1. 

We can now start again from (/i7-. Ay) at time T to extend the solution to [0,2T], 
and so on, to prove the proposition. 

We use the following norm on x M: 

||(/i, A)|| = \\fi\\ + \X\. 
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We note the following estimates: there is a constant C < oo, depending only on (p 
and B such that, for all /i, fi' G and all A, A' G M, 

(2.7) ||L^(/x,A)||<C||(/i,A)ir, 

(2.8) 1|L^(^,A) -L^(^',AO|| < C||(/x,A) - (/i',A')|| (||(^,A)|| + ||(/.', A')||). 

We turn to the iterative scheme. Set {fj,^, A^ ) = (/xq, Aq) for all t and define inductively 
a sequence of continuous maps 

t ^ (/i^", a;') : [0, cx)) ^ Mb X M 

by 

(/xr\ Ar^) = (/io, Ao) + /* L^(/i:, ad ds. 

Jo 

Set 

/„(t) = ii(/ir,Ar)ii 

then fo{t) = /„(0) = ||(/io, Ao)|| < 1 and by the estimate (g^) 



/n+i(t)<l + C / MsYds. 

Hence 



t 

2 







Ut) < (1 - ct)-\ t<c-' 

for all n, so, setting T = (2C)~^, we have 

(2.9) ||(/xr,Ani|<2, t<r. 

Next, set goit) = fo{t) and for n > 1 

g^it) = \\{f,-,X7)-if^r\K~')\l 

By the estimates (|2.8|) and ( ^l9|) , there is a constant C < cxd, depending only on ip 
and -B, such that 

ft 



9n+i{t) <C gn{s) ds, t <T. 
Jo 



Hence, by the usual arguments, (/i". A") converges in x M, uniformly in t < T, to 
the desired local solution, which is also unique. Moreover, for some constant C < oo, 
depending only on ip and B, we have 

\\{f^t,Xt)\\<C, t<T. 

It remains to show that fit ^ for all t. For this we need the following result. 
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Proposition 2.3. Let 

{t,x)^ft{x) : [0,T] X 

be a bounded measurable function, having a bounded partial derivative df/dt. Then 
for all t < T 

I = {df/dt,^,t) + {{ft,0),L^{fit,Xt)). 

Proof. Fix t <T and set [sj„ = {n/t)~^[ns/t\ and [s]„ = {n/t)"^\ns/t]. Then 

{ft, f^t) = {fo, t^o) + [ {df/ds,fiisu)ds+ [ 0),L^(/i„ A,)) rfs 

Jo Jo 

and the proposition follows on letting n — > oo. □ 
For t < T, set 

6t{x) = exp y K{x,y)iJ,s{dy) + Xs(p{x)^ ds 
and define Gt : Ms by 



(/.^^(/i)) = i / UOt){x + y)U+y^BK{x,y)et{x)-%{y)-^^^{dx)^Ji{dy). 

^ J ExE 



We note that Gt{ii) > whenever > and, for some C < cxd, depending only on 
and B, we have 

[[^^(/x)!! < CM\ - < - /i'll (||/i|| + ll/i'll). 

Set /ij = By Proposition ^]3|, for all bounded measurable functions /, we have 



j^{f,h) = {fde/dt,f,,) + ((A,0),L^(/x„A,)) = {f,Gt{jlt)). 
Define inductively a new sequence of measures /t" by setting /i° = /xq and, for n > 

= j2o+ [ Gs{fi'^)ds. 
Jo 

By an argument similar to that used for the original iterative scheme, we can show, 
first, and possibly for a smaller value of T > 0, but with the same dependence, that 
is bounded, uniformly in n, for t <T, and then that — — >■ as — > cxd. 
Since /i" > for all n, we deduce jit ^ ^ and hence /x* > for all t < T. This 



completes the proof of Proposition p.2| □ 



We remark that the arguments used to prove Proposition |2.2| apply with no essential 
change to the case where the coagulation kernel is time-dependent provided that ( p.4|) 
holds uniformly in time. We remark also that, in the iterative scheme 

/i° = /xo, /^r' < /io + / (/is + * /i") ds 

Jo 
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for all n > 0, SO by induction we have 

oo 
k=l 

where /Iq" denotes the n-fold convolution of fiQ. On letting n — ^> oo, we see that, if 
{fit, ^t)t>o is the unique solution to (|2.(j|), then Ht ^ 1q- These remarks will be used 
in the proof of Proposition W^. 



We now fix /xq ^ 3YC^ with /xq > and /io) < oo. For each compact set B C E, 

let 

/i^ = lij/io, = / v5(x)/io (rfx) 

and denote by (;uf , Af )t>o the unique solution to ( p.6|) , starting from (/x^, A^), pro- 
vided by Proposition |2.2| . We shall show in Proposition p.4| that for B C B' we 
have 

/if < fif, {if, /if) + Af > {if, /if) + Af '. 

We shall also show in Proposition p.5| that, for any local solution {vt)t<T of the 
coagulation equation ( |2.1|) , for alH < T 

/if<i/t, (¥.,/if) + Af > (vp,z/i). 



We now show how these facts lead to the proof of Theorem ^TT 
Set /if = Mvub^e IJ't ^iid Xt = limBiE Af . Note that 



{if, fit) = lim(<^,/it ) < (<^,/io) < oo. 



So, by dominated convergence, using ( p.4| ), for all bounded measurable functions /, 

/(x + y)l^+y^BK{x,y)iJ,f {dx)fif{dy) 



and we can pass to the limit in (2^) to obtain 

-rAf, l^t) = - / {f{x + y)- fix) - f{y)}K{x, y)fit {dx)fit{dy) - \t{f^p, fit). 

For any local solution {i't)t<T, for all t < T, 

fJ't<^t, (v5,/it) + Ai > (v?, z^t). 

Hence, if Af = for all t < T, then {fit)t<T is a local solution and moreover is the 
only local solution on [0,T). If {vt)t<T is a strong local solution, then 



{^'^,fis)ds< / {(p"^ , Ug) ds < oo 
Jo 

for all t < T; this allows us to pass to the limit in ( |2.6| ) to obtain 
(2.10) ^Af = Af(v.2,/if) 
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and to deduce from this equation that = for all t < T; it follows that {i't)t<T is 
the only local solution on [0,T). Note that, for any local solution {i't)t<T 

xlx<nJ^t (dx) = / xlxKnJ^O (dx) 
E Je 

+ \ I I {{x + y)lx+y<n- xlx<n-y'^y<n}K{x,y)vs{dx)vs{dy)ds. 
^ Jo Jexe 

Hence, if {i't)t<T is strong and ip{x) > ex for all x, for some e > 0, then by dominated 
convergence the second term on the right tends to as n — > oo, showing that {vt)t<T 
is conservative. 

Suppose now that (v?^,/io) < oo and set T = {ip"^, fj,o)~^. For any compact set 
B C E we have 

TM^^'^) <\l Wix + yr-^ixr-^iyr}Kix,y)f^fidx)f^fidy) < 

^ J ExE 

so, for t < T, 

(¥.^/x,)<lim(y,^/if)<(T-^)-^ 



Hence (|2.10| ) holds and forces A* = for t < T as above, so {fit)t<T is a strong local 



solution. 

If (f"^ is sublinear, then 

(v5^/it) < (v5^/uo) < oo. 

If, on the other hand, K{x,y) < ip{x) + f{y), then 
d 



dt 



{ip\fif) < / ipixMy){ipix) + (^(y))/if (dx)/if (rfy) 



ExE 



< 2(<^,/if)(<^^/if) < 2(<^,/io)(<^',/^n, 

so 

(V5^/i^) < exp{2(v?,/io)t}. 
In either case we can deduce that {fit)t>o is a strong solution. 

Proposition 2.4. Suppose B C B' . Then, for all t > 0, 

/if < /if, {ip,fif) + Af > {ip,fif) + Af . 

Proof. Set 

dt{x)=expj^ (^j^K{x,y)fif{dy) + Xf^ix)^ ds, 
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Note that ttq > and Xo = 0. By Proposition p.3| , for any bounded measurable 
function / 

^(/,vr,) = {f^e/^t,^^f + ((A,0),L^'(/xf , Af) +^''(/if, Af)) 



dt 



fetix + y)K{x,y) 

ExE 

X {U+y(.B'f^f{dx)fif{dy) - U+y(.Bf^f{dx)nf{dy)) 
+ [ fetix)i^ixMy) - Kix,y))fif{dx)ififidy) - fxfidy)) 

J ExE 

+ xt [ fetixMx)fifidx). 

J E 

Also 

:i;Xt = ^ / {ipix) + ip{y)-ip{x + y)}K{x,y){fif'{dx)fif{dy)~fif{dx)fif{dy)). 
^ Jexe 

So (vr^, Xt) satisfies an equation of the form 

^(TTt, Xt) = Ht{iTt, Xt) + {at, 0) 

where Ht : M^' x M — > M^' x M is linear, Ht{TT, x) > whenever (vr, x) > 0, where 
at E M,B' with at > 0, and where we have estimates, for t < 1, 

\\Ht{7i,x)\\<C\\{7i,x)\\ and M < C 

for some constant C < oo depending only on (p and B'. Therefore we can apply the 
same sort of argument that we used for non-negativity to see that tt^ > and Xt ^ 
for alH < 1, and then for alH < cxd, as required. Explicitly, Ht is given by 



((/, a), Ht{7i, x)) = 1 I fOt{x + y)l^+y^BK{x, y) 

^ -'ExE 



X {9t{x)~^7i{dx)nf (dy) + 9tiy) ^nf{dx)7r{dy)) 



+ / fet{x){v{x)^{y) - K{x,y))^if\dx)et{y)-Mdy) 

J ExE 

+ X [ f9t{x)ipix)fxf{dx) 



+ / {v^i^) + ¥^iy) - + y)}K{x, y) 

J ExE 

{et{x)-'7i{dx)fif'{dy) + 9t{y)-'fif{dx)n{dy)) 



2 



and at is given by 



(/> at) = 7: / + y)U+y(LB'\BK{x, y)fif{dx)fif{dy). 

^ Jexe 
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□ 



Proposition 2.5. Suppose that {vt)t<T is a local solution of the coagulation equation 
( |2.1| ), starting from /iQ. Then, for all compact sets B C E and all t <T , 



Proof. Set 



We have to show that tt^ > and Xt > 0. By an obvious modification of Proposition 
2!3| we have, for all bounded measurable functions /, 

|(/,7r,) = {fde/dt,u^-i^f) + {fetlB.L{u,)) - ((/^„0),L^(/xf, Af)). 



By (pT3|), we have 



1 



(V^, ^t) < (V2, J^o) + 7: / + - - ^{y)}K{x, y)us{dx)us{dy) ds. 

Jo Jexe 

On the other hand 

d \ f 

-({if, fif) + Af ) = - / Mx + y)- cfix) - ip{y)}K{x, y) fif (dx) fif (dy) . 
"I ^ Jexe 

So Xt > Pt where 

Pt = 7; / + v{y) - <^{x + y)}K{x,y){us{dx)us{dy) - fif{dx)fif{dy))ds 



2 Jo 



ExE 



Now {nt, Pt) € Mb X M obeys a differential equation of the form 

^(^t^Pt) = Ht{Trt,pt) + (a*, A) 

where Ht : Ms x M — ^ Ms x M is linear, Ht{TT,p) > whenever (7r,p) > 0, where 
ctt > 0, /3t > and we have estimates of the form 

\\Htin,p)\\<C\\{n,p)\\,\\at\\<C,\f3t\<C. 
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It follows that TTt > and > 0, so also ^ as required. Explicitly, Hf is given 
by 

((/, a), Ht{n, p)) = 1 I f9t{x + y)l^+y^BK{x, y) 

^ J ExE 

X {et{x)-^Ti{dx)v^{dy) + et{yY^pf{dx)Ti{dy)) 

+ f A(x)(^(a;)^(y) - K{x,y))et{y)-^uf{dx)7,{dy) 
Jexe 



+ p I fet{x)^{x)v^{dx) 

e 



+ \a j {v?(x) + i^{y) - i^{x + y)}K{x, y) 
^ Iexe 



X {dt{x)-^T^{dx)v^{dy) + et{y)-^fif{dx)7r{dy)) 

and at, Pt are given by 

{f,at)= / fetix){ip{x)ip{y) - K{x,y))ly^Bt^fidx)ut{dy) 
Jexe 

E 



+ {Xt-pt) / f6t{x)ip{x)iy^ (dx) 
Je 

f9t{x + y)U+y^BK{x,y)li^^^y)>^BxBVt{dx)vt{dy), 



ExE 



Pt = l:[ + ip{y) - ip{x + y)}K{x, y)l{x,y)^BxBi^t{dx)vt{dy). □ 

-I ExE 



This concludes the proof of Theorem 2.1 



3. An EXAMPLE OF NON-UNIQUENESS 

We construct in this section an example of Smoluchowski's coagulation equation 
having at least two solutions, both of which are moreover conservative. 
Consider the system of differential equations 

(3.1) ^m„(t) = -A„m„(t)m„4.i(t), ra = l,2, ... 

For any solution we have 



m„(t) = m„(0)exp{-A„ / mn+i{s)ds}. 

Jo 

Assume that m„(0) > for all n, then mn{t) > for all n and t. Note that, if is 
given, and we consider the system restricted to n < A — 1, then m„ is decreasing in 
rriN when A — n is odd, and increasing in m^v when A — n is even. 
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Fix N and consider the case 

m„(0) = 

and 



2-", n^l,...,2N 
0, n = 2iV+l,... 



Xn = 8" for all n. 
Proposition 3.1. We have for all n and t 

m2n{t) > ^m2n(0), 

m2n+i{t) < m2n+i(0)exp{-42"i}. 
Proof. Certainly m2N{t) — m2Ar(0) > |m2iv(0). Suppose that n < N and 

m2n+2{t) > ^m2n+2(0) for all t. 
m2n+l{t) < m2„+i(0) CXp{-A2„+im2n+2(0)t/2} 

= m2n+i(0)exp{-42"t} 

f°° 1 
m2n+i (t) dt < m2n+i (0)4-2" = - 8-2- 

m2n(t) > m2„(0)exp{-8-2"A2n/2} > ^m2n(0). 
Hence the proposition follows by induction. □ 

We denote the solution just considered by m^^ . The same arguments establish 
corresponding inequalities for m?^^^ , where the roles of even and odd are swapped. 
Now we let N ^ oo. For N >n, (t) is decreasing in N and ml^_^_i{t) is increasing 
in N for all n. We set 

m^it) = hm mf (t). 

The integral equation 

W = (0) - Xn f ml'' {s)mlUs) ds 
Jo 

passes to the limit to give 

^nit) = m+{0) -Xn [ m+(s)m++i(s) ds 
Jo 



Then 



so 



so 



so m'^ is differentiable with 



jm+{t) = -A„m+(t)m++i(t). 
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So (m^ : n > 1) solves the original system of equations. The same argument produces 
another solution (m~ : n > 1) given by 

m-{t) = hm m'f+\t). 

N—*oo 

We have m+(0) = m-(0) = 2"" for all n. But 

m^nit) < exp{-42-4}m2-„(0) 

for all n and t, so m"*" 7^ m~. 

We now use the solutions m"*" and irT to construct an example of Smoluchowski's 
coagulation equation having at least two conservative solutions. Let xi, 0:2, . . . be an 
increasing sequence in (0, 00) which is linearly independent over Z. For 

x = Xn^ {hxi + . . . + kn-iXn-i), h, ... , A;.„_i e Z 

we write n{x) = n. Denote by / the set of all such x and define n{x) = if a; ^ /. 
Define K : E x E ^ [0,oo) hy 

j^f^^ y\^(K if {n{x),n{y)} ^{n,n + 1} and n > 1, 
1 otherwise. 

Set 



n=l 



and consider Smoluchowski's coagulation equation 
d ,„ . 1 



dt 



if, y^t)^ \ I {f{x + y) - fix) - fiy)}K{x, y)i,t{dy)i,tidx) 
^ Jexe 



starting from /xq. 

According to the definition made in §2, for a solution, we require 



/ / K{x,y)i^s{dx)i^s{dy)ds < 00 

Jo JBxE 



for all t and all compact sets B C E, and, for all bounded measurable functions / of 
compact support 

(/, /^i) = (/, /^o) + ^ / / {f{x + y) - f{x) - f{y)]K{x,y)iis{dx)^s{dy)ds. 
Consider for n = 1, 2, . . . , 

mn{t) = iit{{x : n{x) = n}). 
Take f{x) — ln{x)=n,x<k s-nd let A; ^ 00 to obtain 

^m„(t) = -A„m„(t)m„+i(t). 
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We deduce that any solution {iit)t>o of the coagulation equation gives rise to a solution 
{mn{t) : n > l)t>o of the system On the other hand, for any solution {mn{t) : 

n > l)i>o of this system, we obtain a solution {fit)t>o of the coagulation equation by 

d X 

— = -(A„_im„_i(t) + A„m„+i(t))/it({x}) + ^>^n-i ^ fit{{y}) fit{{z}) 

y,z&I 

y+z=x 

whenever n{x) = n. Hence the coagulation equation has two distinct solutions {^t)t>o 
and {fit)t>o corresponding to {m^)t>o and (m^)i>o. We now show these solutions 
are conservative. The idea of the proof is to show that the proportion of original 
particles making at least 2n jumps falls off geometrically in n. 

Proposition 3.2. Suppose that has finite total mass density. Then the solution 
{fi^)t>o is conservative. 

Proof. For x = kiXi + . . . + fc„_ia;„„i + x„ G /, define 

{km if m < n 
1 if m = n 
if m > n. 
Note that km : I Z+ is additive. For m < n consider 



Vm,nit) = J km{x)ln{x)<nlJ't{dx). 

Then Vm.n{t) is non- increasing in t. Note that 

km{x)ln{x)<nK{x,y)fit{dx)fit{dy) < 8'^Vm,n{t) < oo 



Ixl 



so, by dominated convergence, since {fit)t>o is a solution. 



^Vm,nit) = -A„m„+i(t) / kmix)ln{x)=ntJ'tidx). 



Note also 



so 



— J km{x)ln(x)=nl^t{dx) < --^Vm,n-l{t) 



km{x)ln{x)=nfJ't{dx) < rm,n-l{t) 

where rm,n{t) = 2^™ — Vm,n{t)- Now for n even 



1 



{t)< (t) / Xrirrin+iis) ds < -rm,n~iit) 

Jo ^ 

so rm n(t) — > as CX3. Hence Vm nit) 2^™ as n — > cxo and hence 
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□ 

We make some remarks on the relation between this example and the results of §2. 
The construction of the example makes it insensitive to the additive structure of E. 
We require very little of the sequence {xn)n>i- By taking x„ ~ 8" we can satisfy the 
condition 

K{x, y) < X 

but we get 

xfio{dx) = oo. 

On the other hand, by taking x„ ~ a", for some a < 2, we get 

x^o{dx) < oo 

but the relation K{x, y) < Cxy for all x,y & E does not hold for any C < oo. Thus, 
however we choose (x„)„>i, we cannot regard {fif)t>o as a strong solution, even in 
small time. This is of course implied also by the uniqueness of strong solutions 
established in §2. 

It would be nice to find an example of this type where the initial mass distibution is 
supported on N. It may be that for integers x„ oo sufficiently fast, the analogous 
equation exhibits the same sort of behaviour. However we have not established 
whether this is true. 



4. HYDRODYNAMIC limit FOR THE STOCHASTIC COALESCENT 



In this section we shall prove some limit theorems for the stochastic coalescent. 
There are two main results. In Theorem iA, generalizing a result of Jeon [Jeo], we 
prove a tightness result for the stochastic coalescent, which implies a general exis- 
tence theorem for solutions of Smoluchowski's equation. Then, in Theorem we 



prove weak convergence of the stochastic coalescent to any strong solution of Smolu- 
chowski's equation. The methods used are mostly standard tools from the theory 
of weak convergence on Skorohod spaces. The problem-specific idea which leads to 



Theorem [4.4| is the construction of a coupled family of particle systems, converging 
to the stochastic coalescent, in direct analogy with the method of §2. A version of 
this idea was also discovered independently by Kurtz ||Kur|| . The case of a discrete 



mass distribution may also be treated using a differential equation approach instead 
of weak convergence: this is simpler and more effective, establishing convergence at 
an exponential rate in the number of particles. The particle system we consider 
has been considered, in various special cases, by many others. In particular, it was 
considered in full generality by Marcus [[Mar68|| and Lushnikov [[Lus78|| . Recall that 
E = (0, oo) and that the coagulation kernel is a symmetric measurable function 
K : E X E [0,oo). 
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Let Xo be a finite, integer-valued measure on E. We can write Xq as a sum of unit 
masses 

m 

1=1 

for some Xi, . . . ,Xm & E. We think of Xq as representing a system of m particles, 
labelled by their masses xi, . . . , Xm- A Markov process {Xt)t>o of finite, integer-valued 
measures on E can be constructed as follows: for each pair i < j, take an independent 
exponential random time Tij of parameter K{xi,Xj) and set T = minj<jTjj; set 
Xj = Xq for t < T and set 

= Xq — £3;, — + ^Xj+xj if T = Tij] 

then begin the construction afresh from Xt- In this process, each pair of particles 
{xi,Xj} coalesces at rate K{xi,Xj) to form a new particle Xi + xj. We call {Xt)t>o a 
stochastic coalescent with coagulation kernel K. 

Denote by d some metric on M which is compatible with the topology of weak 
convergence, that is to say /i) — > if and only if (/, (/, /i) for all bounded 

continuous functions / : — > M. We choose d so that d{^, fi') < — for all 
/i, /i' G M. When the class of functions / is restricted to those of bounded support we 
get a weaker topology, also metrizable, and we denote by do some compatible metric, 
with do < d. 

The following result is a first attempt at proving weak convergence for the sto- 
chastic coalescent. It is less than satisfactory because it does not enable us to show 
uniqueness of limits. We include it here, partly as a warm-up for the more intricate 
arguments used later, and partly because it provides the best result on global exis- 
tence of solutions to Smoluchowski's equation that we know. A version of the result 
where /xq is supported on N and where ip{x) = x has been proved already by Jeon 
Jeo| . 

Theorem 4.1. Let K : E x E ^ [0, 00) be a symmetric continuous function and 
let fiQ be a measure on E. Assume that, for some continuous sublinear function 
(p : E —* {0, 00), 

K{x,y) < (p{x)(p{y), forallx,yeE, 

(fixy^ifiyy^Kix, y) 0, as (x, y) 00. 

Assume also that {ip,fio) < 00. Let {X^)t>o be a sequence of stochastic coalescents, 
with coagulation kernel K . Set X" = n~^X^^i^ and suppose that 

do{LpXQ,Lpno) -> 
as n 00 and that, for some constant A < 00, for all n 

(^,Xo")<A. 
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Then the sequence of laws of on D([0, oo), (M, rfo)) is tight. Moreover, for any 
weak limit point X , almost surely, {Xt)t>Q is a solution of Smoluchowski's coagulation 
equation ( p.l| ) . In particular, this equation has at least one solution. 



Proof. For an integer-valued measure jj, on E, denote by ijl^^'> the integer-valued mea- 
sure on E X E given by 

X A') = fi{A)fi{A') - fi{A n A'). 

(This is simply the counting measure for ordered pairs of masses of distinct particles.) 
Similarly, when ufi is an integer-valued measure, set 

X A') = fx{A)fi{A') - n-^fi{A n A') 

and note that n^/i^") = (njj,)^^\ For a bounded measurable function / on E, set 

= 11 {fix + y)- fix)- fiy)}K{x,y)fx^'^\dx,dy), 
= \l {f{x + y) - f{x) - f{y)YK{x,y)^,^'^\dx,dy). 

ExE 



2 

Then 

M/'" = - - f L^'\x:)u) ds 

Jo 

is a martingale, having previsible increasing process 

(M/-")i= f Q^'\x:){f)ds. 
Jo 

Set M/'" = n'^M-^l" , then we have 



(4.1) = (/,Xo") +n-^M/'"+ f L^-\x:){f)ds, 

Jo 



Since ip is subadditive, we have {ip, X") < A for all n and t. Hence by (|2 

lL(")(x^)(/)|<2||/||A^ 

Q(")(Xr)(/)<4||/fAl 
Assume that |/| < A 1. Then 

K/,^ni<A 
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for all t, so we have compact containment. Moreover, by Doob's L^— inequality, for 
all s < t, 



E 



sup |M/'" - M/'"|' < 4E / g(")(Xn(/)f^^ < 16A'(t - s) 

s<r<t Js 



SO 

E sup \{f,X^-X^)\^<C{{t-sy + n"\t-s)} 

s<r<t 

where C < oo depends only on A. Hence, by a standard tightness criterion, the laws 
of the sequence (/, X") on D([0, oo),]R) are tight. See for example [[EK86|] , Corollary 
7.4. We note the bound 

ii(^Ai)xrii<(vp,xr)<A 

for all t. Hence we can apply Jakubowski's criterion ||Jak86|| to see that the laws of the 
sequence (y? A on -D([0, oo), M[o,oo]) are tight. By consideration of subsequences 
and a theorem of Skorohod, see, for example. Pollard [Pol84|| , Chapter IV, it suffices 
from this point on to consider the case where {(f A converges almost surely in 

D{[0, oo), M[o,oo]), with limit (99 A 1)X say. We denote also by X the process in Me 
obtained by restriction of measures. Note that 

\\X^-Xl\\<S/n 

so X G C([0, 00), M^;). Moreover ip^ = ipl{o,s] is subadditive, so 

xn < {ip', X") < {if', /io) + 1 - /.o) i 

and so, given e > 0, we can find 5 > so that 

sup sup(v9'',X") < e. 

n t 

Given a continuous bounded function / : ii^ ^ M of bounded support, we can write 
/ = /i + /2 where fi is continuous of compact support and /2 is supported in (0, 6) 
with II/2II < ll/ll. Then 

limsup snp{ipf,X: - X,) < lim snp{ipfi,X: - X,) + ||/|| sup snp{ip'x:) < e 

rt— >oo s<.t n~KX) ^<;^ ^ g 

Since / and e were arbitrary, this shows that 
(4.2) s\ipdo{cpX^,ipXs) -^0 a.s. 

s<t 
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We now wish to pass to the hmit in ( [4.1| ). Let us suppose for now that / : E 
is continuous and of compact support B. Then, as n oo 

E(sup |n-5Mf'"n < lE(M^'")t < ^ 0, 

s<t n 11 



l(z,_z,W)(x;)(/)| = i 



{f{2x)-2f{x)]K{x,x)X:{dx) 



< ^ / ^{xfx:{dx) 

IBU2B 



n 



n 

Hence it will suffice to show that as n ^ oo 

(4.3) sup|(/,L(X;)-L(X,))| ^0 a.s. 

s<t 

where we recall that 

(/, L{^)) = \ j {f{x + y)- fix) - f{y)}K{x, y)fx{dx)fx{dy). 

Given 5 > and N < oo, we can write K = Ki + K2, where Ki is continuous of 
compact support and where < K2 < K and K2 is supported on 

Fi U F2 U F2 = {{x,y) ■.x<5}U {ix,y) : y < 5} U : > A^}. 

Then, with an obvious notation, by 



sup|(/,Li(X;)-Li(X,))|^0 a.s. 

s<t 

whereas, for K2, we use the estimates 

WKlp^fi® = WKlp^fJ' ^ fJ-W < (v2,yu)(v5'^,/i), 
Iliriirg/iO/ill < /5Jv(v^,/i)^ 

where (3n = sup^^^^y^^^^ (p{x)-^ip{y)~^K{x,y). Now, 

and, given e > 0, we can find 6 and so that 

{^\X^)<leA~\ for all n. 

Then 

\{f,L2{X:))\<s, \{f,L2{X,))\<e 
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for all n and t. Hence 

limsup sup |(/, L{X:) - L(X,))| < 2e. 

n— >oo s<t 

But e was arbitrary, so ( [4.3D is proved. Hence we can let n — > oo in ( ^.11 ) to obtain 

(/,X,) = (/,Xo)+ [\f,L{X,))ds 
Jo 

for all continuous functions f : E ^ M. of compact support. By using the bounds 
(|2.4| ) and {(p,Xt) < A, and a straightforward limit argument, we can extend this 
equation to all bounded measurable functions /. In particular, almost surely, X is a 
solution of Smoluchowski's equation, in the sense of §2. □ 



A corollary of Theorem of |4.1| is that, whenever we know Smoluchowski's equation 
has at most one solution, then, under the hypotheses of Theorem O, we can deduce, 
for all t 

sup(io(X",/x^) 



in probability as n — >■ oo, for the solution {fJ,t)t>o provided by Theorem [4.1| . However, 
we can only prove uniqueness of solutions in the presence of a strong solution. So we 
prefer to formulate our main limit result. Theorem [4.4| , in that context, when a new 
approach allows certain other hypotheses to be relaxed. 

For the remainder of this section we will assume that we have chosen a continuous 
sublinear function ip : E ^ {0, oo) and that K satisfies 

(4.4) K{x,y) < ip{x)ip{y) for all x, ?/ G -E. 

Our further analysis of the stochastic coalescent will rest on an approximation by 
a coupled family of Markov processes {X^ , A^)t>o, indexed by sets B C E which 
we now describe. Each process {X^)t>Q will take values in the finite integer-valued 
measures on E, whereas (Af)t>o will be a non- decreasing process in [0, oo). Let us 
suppose given initial values [Xq, Aq), for all B, such that X^ is supported in B and 
such that, whenever B (1 B' 

Xi < Xf , Xi) + Ao^ > Xf ) + Af . 

Write Xq = Xq as a sum of unit masses 

m 
i=l 

where xi, . . . ,Xm^E. There is a unique increasing map 

B ^ I{B) C {1,... ,m} 



such that 



i€l(B) 



SMOLUCHOWSKI'S COAGULATION EQUATION 23 

Set 



Note that decreases as B increases and that = Aq > 0. For i < j, take 
independent exponential random variables Tij of parameter K{xi,Xj). Set Tji = Tij. 
Also, for i ^ j, take independent exponential random variables Sij of parameter 
ip{xi)Lp{xj) — K(xi, Xj). We can construct, independently for each i, a family of inde- 
pendent exponential random variables S^, increasing in B, with Sf having parameter 
ip{xi)i'^ . Set 

and note that is an exponential random variable of parameter 



For each B, the random variables 

{T,„Tf eIiB),t<j) 

form an independent family, whereas, for i G I{B) and j ^ I{B), we have 

< T- ■ 
-'-I — -'jj 

and for B C B' and all i we have 

rriB ^ rjiB' 
i — i ■ 

Now set 

r= (mmTij) A (minT^*^). 

i<j i 

We set (Xf , Af ) = {X^ , A^) for t < T and set 

{{X^ - e^. - e^. + Af ) if T = T^^-, i, j e /(fi), + x^- G 

(Xg^ - e^. - So;,-, + ip{xi + a;^)) if T = Tj^, i, j G I{B),Xi + Xj ^ B 

W - Af + ^(o;,)) if T = i;^, z G I{B) 

(X^,Af) otherwise. 

It is straightforward to check that X^ is supported on B and, for B C B' 

X^ < X^', {if, X^) + A^ > (if, X^') + Af . 

Hence we can repeat the above construction independently from time T, to obtain 
a family of Markov processes {X^ , Af)t>Q such that Xf is supported on B and for 
B CB' and all t 

(4.5) Xf < Xf, Xf ) + Af > {cp, Xf) + Af '. 
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At the outset, we assumed that both and A^ were given, for all B. Prom now 
on we shall suppose simply that Xq — X^ is given and take 



Of course these relations do not remain valid as time evolves. 

For each fixed S, the process {X^ ^ A^)t>o may be regarded as a finite state-space 
Markov chain having three sorts of transition. Each pair of particles Xi,Xj in X^ is, 
at rate K{xi, xj), removed; if Xi + xj e B, the merged particle is added to Xf, if not, 
(p{xi + Xj) is added to Af . Also, each particle Xi in X^ is, at rate ip{xi)A^, removed 
and ip{xi) added to Af . In particular, for the choice of initial values made above, 
Af = for all t and Xt = Xf is simply the stochastic coalescent with coagulation 
kernel K with which we began. 

We now proceed to identify some martingales associated with (Xf, Af )t>o. Recall 
that, when nfj, is an integer- valued measure on we denote by //^"^ the measure on 
ExE characterized by 



Given an integer- valued measure ji on E and given A > 0, define for any bounded 
measurable function f on E and for a e R 



^0 — ^bXq, 




{cplBc,Xo). 



^"\AxA')^l^{A)i^{A') 



n-'fx{Ar)A'). 



L^'«(/.,A)(/,a) = ((/,a),L^'W(/.,A)) 




f{x)-f{y)}K{x,y)f,^'\dx,dy) 




and 



Q^'W(/.,A)(/,a) 




f{x)-fiy)}'Kix,y)i,^'\dx,dy) 




Then, for all / and a, 



(4.6) 



Mt = (/, ) + aAf - (/, Xi) - aA, 




is a martingale, having previsible increasing process 



(4.7) 
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Recall from §2 that, for B C E compact, we denote by the space of finite signed 
measures supported on B and we define iM^xR— ^M^xMby the requirement 

((/,a),L^(/i,A)) = 

J / {f{x + y)l^+y(.B-a'^{x + y)l^+y^B-f{x)-f{y)}K{x,y)fi{dx)fi{dy) 
^ Jexe 

+ A / {aip{x) — f{x)}(p{x)iJ,{dx) 



for all / and a. 

Fix a measure /io on E with {(p,fio) < oo. Set 

/i^ = Is/iO, = (v3lBc,/io). 

Recall from §2 that, for each compact set B, the equation 



(4.8) (/if, Af ) = (/.o", A^) + / L-(/.f , Af ) 

Jo 

has a unique solution, which is a continuous map 

t ^ (/if, Af ) : [0, oo) ^ X M+. 

Consider now a sequence of integer- valued measures X^. For each n, denote by 
(X")t>o and (X^ '")t>o the stochastic coalescent and the coupled family of ap- 
proximations constructed above, starting from Xq. Set 

Xf = 

{X^'\Afn=n-\X^ll,AX)- 

We shall need a mild continuity condition on K. Denote by S{K) (1 E x E the 
set of discontinuity points of K and by /i™ the ra-th convolution power of fiQ. Our 
assumption is that 

(4.9) {fi*o''f\S{K)) = 0, for all n > 1. 

This condition is verified, in particular, when S{K) has Lebesgue measure zero and 
/io is absolutely continuous. 

For the purposes of the next proposition, we also need an analogous condition on 
the compact set B: 

(4.10) I^Ti9B) = 0, for all n > 1. 

This condition is verified, for any given /io, for all but countably many closed intervals 
in E. 

Proposition 4.2. Assume conditions (p. 41), (|2.5| ), ( ^4.9|) , ([4.10|) . Suppose that 

d(Xo^'",/i^)^0, |Af'"-Ao^|^0 



26 JAMES R. NORRIS 

as n —>■ oo. Then, for all t > 

sup d(Xf /xf ) ^ 0, sup |Af - Af I ^ 

s<t s<t 

in probability as n —>■ oo. 

Proof. Set A = sup„(v9,Xq) and note that A < cxd. By rescaling (f4.6|) and ([4.7|) , we 
see that, for all B, all bounded measurable functions / and all a G M 

(4.11) 

Mr = Xf + aAf - (/, Xo^'") - aA^'" - L^'(")(Xf Af '")(/, a) t/.) 

is a martingale, having previsible increasing process 

(4.12) (M")i= rQ^'(")(Xf",Af'")(/,a)ds 

^0 

where 

L^'(")(;,,A) = n"2L^'«(n^,nA), 

g^'W(/i,A) = n-2g^.(i)(n^,riA). 
There is a constant C < oo, depending only on B,A and such that 

\L^iX^'-,Afnif.a)\<Ci\\f\\ + \a\), 

_ L^.("))(xf'^ Af '")(/, a)| < Cn~\\\f\\ + |a|), 

|Q^'(")(Xf'", Af '")(/, a)|<C(||/|| + |a|)^ 

Hence by the same argument as in Theorem [4.1|, the laws of the sequence (X", A") are 
tight in -D([0, cx)), M^xM). Indeed, similarly, the laws of the sequence (X", A", J", J") 
are tight in D{[0,oo),Mb x M x Mbxb x ^bxb), where 

indx,dy) = K{x,y)l,+y^BXl'{dx)X^{dy), 

JJ}{dx,dy) = K{x,y)l^+y^BXJ^idx)XJ'{dy). 

Denote by (X, A, /, J) some weak limit point of this sequence, which, by passing to 
a subsequence and the usual argument of Skorohod, we may regard as a pointwise 
limit in D{[0, oo),Mb x R x M^bxE x ^Bxb)- Then there exist boimded measurable 
functions 

/, J : ^] X [0, oo) X X ^ [0, oo) 
symmetric on B x B, such that 

It{dx,dy) = I{t,x,y)Xt{dx)Xt{dy), 
Jt{dx,dy) = J{t,x,y)Xt{dx)Xt{dy) 

in Mbxb and such that 

lit^x^y) = K{x,y)U+y(.B 
J{t,x,y) = K{x,y)U+y^B- 
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whenever [x, y) ^ S{K) and x + y ^ OB. Moreover we can pass to the hmit in ( |4.11|) 
to obtain, for all continuous functions / and all a G M, for alH > 0, almost surely 

((/,a),(Xi,Ai)) = ((/,a),(Xo,Ao)) 



1 



+ 7j/ / {f{x + y)~f{x)-f{y)}I{s,x,y)Xs{dx)XMy)ds 
^ Jo Jexe 



(4.13) 1 



+ + y) - f{x) - f{y)}J{s,x,y)Xs{dx)Xs{dy)ds 

^ Jo J E>(E 



+ As {a(p{x) - f{x)}ip{x)Xs{dx)ds 



By the remarks following the proof of Proposition p.2| , this equation forces Xt ® Xt 
to be absolutely continuous with respect to 



oo 



n=l 



Hence by the assumptions ([4.9|) , (|4.1CI| ), we can replace I(t,x,y) by K{x,y)lx+y(zB 
and J{t,x,y) by K{x,y)lx+y^B in ( |4.13|) . But this is now equation ( [4 .81) which has 
a unique solution (/if, Af )f>o. We have shown that the unique weak limit point 
of (X",A"') in D([0, oo),Mb) x M is the continuous deterministic path (/if,Af)t>o, 
which proves the proposition. □ 

We consider now the special case where /iq is a probability measure on N = 
{1,2,...}. Here we can replace the method of weak convergence in Proposition 
4.2 by a more direct approach using differential equations. The benefit in this ap- 
proach, besides greater transparency, is that we can establish a rate of convergence, 
which is in principle computable. This would be needed if one wished, in practice, 
to assess whether Smoluchowski's equation provided a tolerable approximation to 
the stochastic coalescent. Since the stochastic coalescent already makes a mean- 
field approximation — it is assumed we neglect spatial variations in the particle mass 
distribution — we are effectively assuming there is some external spatial mixing. The 
relevant particle number n is then the number of particles in the largest region which 
is mixed to equilibrium in unit time. 



Proposition 4.3. Assume conditions ( p.4| ) and ( |2.5| ). Suppose that /iq is supported 
on N. Let B be a finite subset ofN. Then there is a constant C < oo, depending only 
on K , yj, hq, and B such that, for all n>l, for X^ = n~^X^^i^ and Af = n~^A^^i^, 
for aU0<5<t, %f 5o = \\X^ - + |A^ - < 1 then 

P (sup{||Xf - /.f II + |Af - Af 1} > (5o + 5)e^*) < Ce-^^^^^*. 
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Proof. We regard {X^ , A^)t>Q as taking values in the finite-dimensional vector space 
X M. Recall from the proof of Proposition |4.2| that 

(4.14) M, = v/^{(Xf ,Af) - (Xo'',Ao^) - T L^'(")(Xf , Af ) rf.} 

Jo 

is a martingale, having previsible increasing process 

Jo 

We recall a form of the exponential martingale inequality for martingales M whose 
jumps are bounded uniformly by A G [0, oo) and which have a continuous previsible 
increasing process (M): for all 6* > we have 

P(supMi > 6 and {M)^ < < exp{-0(5 + -Oh^^e}. 
t 2 

Let f : B ^ M. and a G M be given, with ||/|| + \a\ < 1. Consider the martingale 

M/''^ = ((/,a),M,). 

The jumps of M-^'" are bounded uniformly by Cn^^ for some C < oo, depending on 
ip and B. The process {M^'"-) is continuous and satisfies 

(M-^'")i < Ct 

for some C < oo, depending on ip, B and fiQ. So, by the exponential martingale 
inequality 

P(supM/''^ >6) < exp{-96 + -COHe^'^^^}. 

s<t 2 

Assume that 5 < and take 9 = 5/{3Ct). Then CO/y/n < 1/3, so e^^/^ < 3/2 
and so 

s<t 

Since B is finite, we deduce that for some C < oo, depending on ip, B and fiQ 

P(sup||M,|| >S) < Ce-^'/^* 

s<t 

whenever 6 < -Jnt. Hence 

P(sup \n-^Ms\ >b)< Ce""-^'/^* 

s<t 

whenever 5 <t. Note also the estimate 

_ , Af)|| < Cn-\ for all t > 0, 
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for some C < oo, depending on ip and B. Set Yt = {X^ , kf) — (/if, Af ) and subtract 
the equations ( [4.14| ), ( |2.6| ) to obtain 

Yt = Rt+ [ Lf{Y.,)ds 
Jo 

Rt = Yo + n-'^Mt+ - L^'("))(Xf , Af ) ds 
Jo 



where 



and where 

((/, a), Lfifi, A)) = i / {f{x + y)lx+yeB + a^{x + y)lx+y^B - f{x) - f{y)} 

Jexe 

xK{x,y){Xf + fif){dx)fi{dy) 
+ (Af + Af ) / (a^ix) - f{x)Mx)fi{dx) 

J E 

+ A / (a^(x)-/(a;)Ma;)(Xf + /xf)(cia;). 
Je 

Note the estimate 

||Lf(/^,A)||<C||(/^,A)||/2 

where C < oo depends on B and /io- Set g{t) = supj,<j and r{t) = supj,<( ||-Rs||- 
Then 

9{t)<r{t) + ^-C J\{s)ds 

so g{t) < r(t)e*"*/^. Now, for 5 < t, we have 

P(r(t) > ^(0) + 6/2 + Ct/n) < Ce"'"'''/^*. 
We may assume that C > 4. If 5^ < Ct/n we have nothing to prove. Otherwise 

Ct/n <S^< i<5e^*/2 

so 

Hgit) > {6o + 6)e^') < P(r(t) > ^(0) + 6/2 + Ct/n) < Ce""'''/^*. □ 
Here is the main result of this section. 

Theorem 4.4. Let K : E x E —>■ [0, oo) be a symmetric measurable function and let 
fiQ be a measure on E. Assume that 

(/iS")®'(5(ir)) = for all n> I 

where S{K) denotes the discontinuity set of K . Assume also that, for some contin- 
uous sublinear function ip : E {0, oo), 

K{x, y) < ip{x)ip{y) for all x,y E E 
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and that (v5,yUo) < oo and ((y9^,/io) < oo. Denote by {fit)t<T the maximal strong 
solution provided by Theorem ^. i| . Let {X^)t>o be a sequence of stochastic coalescents, 
with coagulation kernel K . Set X" = n~^X^_i^ and suppose that 

d{<^X^, <^fio) 

as n oo. Then, for all t <T 

sup d{ipX^, ipfis) 

s<t 

in probability, as n oo. 

Moreover, if /ig is supported on N, then for all t < T and all 6 > there are 
constants Sq > and C < oo, depending only on K, hq, ip, t and 6, such that, for all 
n, 

||v^(Xo"-/io)|| <5o 

implies 

P(sup||<^(X,"-/i,)|| >5) <e-"/^. 

s<t 

Proof. Fix (5 > and t < T. Since {fit)t<T is strong, we can find a compact set B 
satisfying ( [4.10D and sucli tliat Af < 5/2. Now 

so 

d{X^^-,f^^)-^0, |Ao^'"-A^|-.0. 

Hence, by Proposition 



sup d(Xf /if) 0, sup |Af - Af I ^ 

s<t s<t 

in probability as n oo. Since {/if : s < t} is compact, we also have 



suprf(y.Xf^y./if)^0 

s<t 



in probability as n — oo. By Proposition and by ( [4.5|) , for s < t 
M^^s - /if) II = /i. - /if) < Af < Af < 6/2 

Mx: - xf '")ii = {^,x: - xf") < Af'" < Af'" 

<Af + |Af'"-Af| 

B,n \B| 



< S/2 + |Af - A 



Now 

d{^x:, ^/i,) < iiv.(x: - xf '") II + rf(^xf ^/if ) + ii^(/i, - /if) 

<5 + d(^Xf'",^/if) + |Af'"-Af| 
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SO 

as oo, as required. 

In the case where n is supported on N, we can argue similarly, replacing the weak 
metric d by the total variation norm and replacing Proposition by Proposition 
O, to arrive at the desired conclusion. □ 



Corollary 4.5. Let K , /io, and Lp he as in Theorem \4.4 - Assume in addition that /xq 
is a probability measure and that Xq is the empirical distribution of a sample of size 
n from fiQ. Then, for all t <T 

sup d{ipX^, v9/i^) -> 

s<t 

in probability, as n ^ oo. Moreover, if /io is supported on N, and if {e°"^, fio) < oo for 
some a > 0, then, for allt <T and all 5 > 0, there is a constant C < oo, depending 
only on K , /io, ^, t and 5, such that for all n 

P(sup||^(Xr-/i,)||>5)<e-"/^. 

s<t 

Proof. For general /io, it suffices to note that 

rf((^Xo", <^/io) 

almost surely as n — > cxd, by the strong law of large numbers, and to apply Theorem 



Suppose now that /io is supported on N. We have 

||(^(Xo"-/io)|| < 2((^l(jv,oo),/io) + |(<^l(jv,oo),^o"-/io)| + (<^l(o,iv], |^o-/"o|) = h+h+h- 

We can choose so that Ji < 5q/2. Then by standard exponential estimates we can 
find C < oo, depending on fio,ip,N and 6o, such that 

On combining this estimate with that found in Theorem E^, we deduce 



P(sup||¥.(X;-/i,)||><5)<2e-"/^ 



s<t 

as required. □ 
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